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The microscopic nature of the hidden order state in URu2Si2 is dependent on the low-energy 
configurations of the uranium ions, and there is currently no consensus on whether it is predomi- 
nantly 5/^ or 5/"^. Here we show that measurement of the basal-plane nonlinear susceptibility can 
resolve this issue; its sign at low-temperatures is a distinguishing factor. We calculate the linear 
and nonlinear susceptibilities for specific 5/^ and 5/^ crystal-field schemes that are consistent with 
current experiment. Because of its dual magnetic and orbital character, a magnetic non-Kramers 
doublet ground-state of the U ion can be identified by Xi{T) oc X3"(^) where we have determined 
the constant of proportionality for URu2Si2 . 



I. MOTIVATION 

Despite tremendous experimental and theoretical ef- 
forts over the last quarter century, the nature of the or- 
dering ("hidden order") at Tq ~ 11. in the actinide 
heavy fermion URu2Si2 remains unresolved [l]. Clear 
signatures of Fermi liquid behavior [TJ [2] above Tq and 
the sharp anomalies in thermodynamic properties [T]-[4] 
at To suggest a Fermi surface instability and thus itiner- 
ant behavior [5 . However, the large entropy of conden- 
sation (-^ ~ 0.3i?ln2), the Curie- Weiss susceptibility at 
room temperature [l] |2] and the strong Ising anisotropy 
observed in the magnetic responses [31 2] all point to a 
localized origin of the key magnetic degrees of freedom 
that contribute towards the formation of the hidden or- 
der state. 

A key outstanding issue for URu2Si2 is to determine 
whether its U ion ground-state configuration is predom- 
inantly 5/^ or 5/^. Specific heat, high-temperature 
susceptibility and even photoemission measurements 
are unable to resolve this issue Inelastic neutron 
scattering experiments [6l [7] favor 5/^ whereas EELS 
measurements [8 are consistent with 5/^. While tech- 
niques like inelastic x-ray and neutron scattering can suc- 
cessfully resolve the crystal field configurations of 4/ ions 
like cerium and praseodymium ilOj, these techniques 
fail for 5/ ions like uranium due to significantly broader 
crystal field levels. Resolution of the ground state con- 
figuration of the uranium ion in URu2Si2 would provide 
an important constraint on the theoretical description of 
its hidden order. 

Here we present the nonlinear susceptibility (xs) as 
a probe of the nominal valence of the uranium ion in 
URu2Si2. In cubic UBeis, X3 provided strong evidence 
that the low-lying magnetic excitations are predomi- 
nantly dipolar (5/^) rather than quadrupolar (5/^) in 
character [1 2 . We generalize this technique to the tetrag- 
onal URu2Si2 system, where the sign of X3 distinguishes 
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Figure 1: Three possible crystal- field schemes for URu2Si2- 
(A) The 5/^ Fs non-Kramers doublet ground-state, with all 
excited crystal- field states degenerate at 1200 K. (B) Two 
finely-tuned low-lying singlets scenario (with the two dou- 
blets found from DMFT[TT]), with all other states degenerate 
at 800 K to reproduce the Van Vleck susceptibility. (C) The 
5/^ Kramers doublet ground state, with a nearly Ising doublet 
ground state, Fe = -.91| ±7/2) + .42| T 1/2) + .05| t9/2) cho- 
sen to best represent the measured Ising susceptibility, with 
another, excited Fe doublet at 400K and all other excited 
doublets above lOOOK. 



between the two nominal valences 5/^ and 5/^. The lin- 
ear and the nonlinear susceptibilities are calculated for 
specific 5/^ and 5/^ crystal-field schemes consistent with 
recent experiment. In addition we discuss the dual nature 
of the magnetic non-Kramers doublet, F5, that is mag- 
netic along the c-axis but quadrupolar in the basal plane. 
As both magnetic (xi) and quadrupolar (xs") suscepti- 
bilities originate from splitting the same doublet, they 
will have the same temperature-dependence even outside 
the single-ion regime. 
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II. TIME-REVERSAL PROPERTIES OF 5f AND 
5/^ CONFIGURATIONS 

The two valence configurations 5/^ and 5/^ are clearly 
distinguished by their properties under time-reversal. 
5/^ is always a Kramers doublet; 5/^ can either involve 
low- lying singlets or, in the tetragonal environment ap- 
propriate for URu2Si2, it can form a Fs doublet formed 
from two degenerate orbitals \x) and \y) that are crystal 
rotations of one another. In a tetragonal environment, 
these two states are degenerate, so one can construct a 
magnetic basis |±) = \x) ± Such a "non-Kramers 

doublet" possesses a dual orbital and magnetic charac- 
ter; and it can be distinguished from a Kramers dou- 
blet through its time-reversal properties. A non-Kramers 
doublet time reverses according to 

^i±) = i±), (1) 

whereas a Kramers doublet transforms with an additional 
and all-important minus sign 

^|±) = ±|T). (2) 

Consider next the case of two almost-degenerate low- 
lying singlets \a) and \b) that are not related by sym- 
metry, separated by a small crystal field splitting A, 
as in the scenarios considered in [11 or \TT. At tem- 
peratures and fields 5, T » A, one can treat this 
almost-degenerate set of singlets as an accidentally de- 
generate non-Kramers doublet, defined by the magnetic 
states |±) = \a) ^ i\b). In this way, the properties of 
a 5/^ and 5/^ configuration at temperatures and fields 
large compared to their splitting can be considered as 
Kramers or non-Kramers doublets respectively. 

A distinctive property of non-Kramers doublets is their 
Ising symmetry, a consequence of their time-reversal sym- 
metry. The matrix elements of a time-reversed operator 
are given by its complex conjugate [M] 

{a\A\l3y = = {a\eAe-^\P), (3) 

where Q\a) = \a) is the time-reversal of state \a). If 
we apply this relation to the off-diagonal matrix element 
{-|J_|-h), then 

(-|J_|+)* = (+|jl|-) = (-|6J_e-i|+). (4) 

Now on the left-hand side, since we obtain 

(— |J_|+)* = (+|J+|— ). The action of time-reversal on 
J_ = — ijy reverses the signs of the angular momen- 
tum operators, and as an anti-unitary operator, takes 
the complex conjugate of the coefficient — z; it thus fol- 
lows that OJ_0~^ = — J+. The time-reversal of the 
non-Kramers states is given by |— ) = |+) and |+) = | — ) 
so that on the right-hand side of Q , (-|6J_B"^|+) = 
— (+1 J+|— ). Comparing both sides, it follows that 

(— |J_|+) = — (+|J_|— ) =0, (non-Kramers doublet) 

(5) 



with a corresponding relation for J+. It follows that 
the off-diagonal matrix elements of J± vanish for a 
non-Kramers doublet, giving rise to a resulting Ising 
anisotropy that is independent of crystal structure. 

By contrast, for a Kramers doublet, the additional 
minus sign in the transformation of the states removes 
the constraint on the off-diagonal matrix elements. To 
see how this works for URu2Si2, suppose that the f- 
configuration is a 5f^ Kramers doublet predominantly in 
an I ±7/2) state. The presence of a tetragonal symmetry 
will add and subtract units of ±4fi, so that the crystal 
field ground-state will take the form 

|±) = a|±7/2)+6|Tl/2>+c|T9/2) (6) 

In this case, |(— |J+|+)p = 56^ + 6ac, so perfect Ising 
anisotropy is only obtained when -\-6ac = 0, i.e if the 
tetragonal crystal fields that mix the two configurations 
are fine-tuned to zero. 

However if the ground-state is a non-Kramers doublet 
of the form 

|±) = a|±3)+6|Tl> (7) 

then an Ising anisotropy follows for arbitrary mixing be- 
tween the I ± 3) state and the | ^ 1) state. Thus the 
observed Ising anisotropy in URu2Si2 either results from 
a finely-tuned 5/^ state, or from a 5/^ state with a real 
or an effective non-Kramers doublet ground-state. The 
challenge is to distinguish these scenarios. We now show 
that this can be done by measuring the nonlinear suscep- 
tibility to a transverse field (perpendicular to the Ising 
axis) . 

Kramers and non-Kramers states differ in their re- 
sponse to a transverse field. If we integrate out the 
high-lying crystal field excitations in the presence of a 
transverse field in the x-direction, the remaining effec- 
tive Hamiltonian will contain off-diagonal terms of the 
form 

Now for a Kramers doublet, applying we obtain 

= -{+\neff{-B)\-) (9) 

so the off-diagonal matrix elements of l-Leff{B) must be 
an odd function of 5, taking the form 

AEK{B)=gB^^jB\ (10) 

Here 7 ~ 1/Ayy is a consequence of third order pertur- 
bation theory, where Ayv is the gap to excited crystal 
field states, and 7 can have either sign. Physically, this 
means that a Kramers doublet can develop dipole and 
octupole components, but has no quadrupolar response. 
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By contrast, for a non-Kramers doublet, applying ([3|, we 
obtain 

{+\neff{B)\-) = {+\eHeff{B)&-'\-) = {+\neff{-B)\- 

(11) 

so the off-diagonal matrix elements of 1-L are even in field, 
where the leading order term is quadrupolar 

AEnk{B) = ^qB^ (12) 

Here q ~ 1/A, since this results from second-order per- 
turbation theory. 

The nonlinear susceptibility, xs^ is defined as the cubic 
term in the magnetization 

M = xi5 + (1/3!)X35^ + .... (13) 

in the direction of the applied field (5), so that xs = 
d^M/dB^ = -d^F/dB^. If we take the the free en- 
ergy F{B) = -Tln^2cosh[/3A£;(5)] 
temperature expansion, we find 

ij^ Aj?(E) y 

^ ^ I I /\ hA M ) I 

(14) 



F{B) = -Tln|2 

1 [mB)y 



and do a high- 
4! V T J 



1 [AE{B)y 
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Using expressions (10) and 

NK _ 3(7^ 

Xs 



T3 
we obtain 



(Non-Kramers doublet) 
(Kramers doublet). (15) 
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We note that the quadrupolar response of the non- 
Kramers doublet leads to a positive Xs^ > 0^ whereas its 
Kramers counterpart is negative at the lowest tempera- 
tures where the 1/T^ term dominates. Thus the different 
time-reversal properties of a 5/^ Kramers doublet and a 
5/^ non-Kramers doublet (or an effective non-Kramers 
doublet composed of two low-lying singlets) can be ex- 
perimentally distinguished by the sign of of the nonlin- 
ear susceptibility at low temperatures. However, it is 
important that the U remain in the single-ion regime. 
If 7 is positive there is a crossover temperature scale, 

Tx ~ ^/If ' where Xs changes sign (but if 7 is negative 

then X3" < for all T). Therefore, provided the single- 
ion regime continues to sufficiently low temperatures, one 
can identify the magnetic state of the U ion from the sign 
of the basal-plane nonlinear susceptibility. This experi- 
ment may require dilution studies to extend the single-ion 
regime down to sufficiently low temperatures. 



III. X3 FOR SPECIFIC 5f AND 5/^ SINGLE-ION 
SCHEMES 

In order to illustrate the above results in detail, we 
have calculated the linear and nonlinear susceptibili- 



ties within three distinct crystal- field schemes: (i) non- 
Kramers doublet, (ii) a "pseudo" -non-Kramers doublet 
given by two closely spaced singlets and (iii) a Kramers 
doublet. These three scenarios are chosen to reproduce 
the high temperature behavior: the Ising anisotropy, mo- 
ment and van Vleck susceptibilities. We assume that the 
maximum in Xi[2^ results from the Kondo effect, and 
thus we are not considering "three-singlet" crystal field 
schemes [13]. 

The three crystal-field scenarios that we consider here 
(cf. Fig. 1) are: 

• (A) Non-Kramers doublet: (J=4) T^: 



|r5±)=a|±3)+6|Tl), 



(16) 



where a, b are adjustable parameters. For simplic- 
ity, all other crystal field states (five singlets and 
one doublet) have been placed at 1200K to repro- 
duce the magnitude of the van Vleck susceptibility; 
however, the behavior of the linear and non-linear 
susceptibilities are relatively insensitive to the ex- 
act arrangement of the excited states. 

• (B) Finely Tuned Singlets : 5/^ (J = 4) with two 
low- lying singlets separated by a gap, A = 36K. 
We use the singlets, r2 and Fi proposed by Haule 
and Kotliar^llj: 



|r2) = 
iri) = 



'1 
2 

cos ( 



V2 



(|4)- 
-(|4> 



4)) 

-4))-sin</.|0), 



(17) 



where (j) = .237r. Again, all excited crystal field 
states have been placed at the same energy, here 
800K, to reproduce the Van Vleck term, and again 
the results are quite insensitive to the exact ar- 
rangement of the excited states. 

• (C) Kramers doublet: bf (J = 9/2) F^^^ 

|F^'^±) =a\± 7/2) + 6| T 1/2) + c| T 9/2), (18) 

where b and c describe the small admixture of 
I ^ 1/2) and | ^ 9/2) chosen to reproduce the 
anisotropy of the susceptibility at room tempera- 
ture. We require that we simultaneously satisfy 
the requirement of perfect Ising anisotropy and the 
correct c-axis moment of 3.5/1^; as we discuss in the 
Appendix, for the Ising state with 56^ + 6ac = 0, 
this is only satisfied with a low-lying crystal-field 
excitation at 110 K that is inconsistent with cur- 
rent neutron scattering results p!^ . We have chosen 
a compromise solution to illustrate the 5f^ config- 
uration (see Appendix for details) with a crystal- 
field scheme that is compatible with experiment [1. 
This doublet is not a pure Ising doublet, but has 
a transverse moment |(— |J+|+)| = .55/iB, giving a 
c-axis anisotropy of 1/50. 
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Based on the high-temperature hnear susceptibihties 
(Fig. |2|, these three scenarios are practicahy indis- 
tinguishable. The nonhnear susceptibihty (Fig 3), Xa"^ 
clearly distinguishes the three cases: the non-Kramers 
doublet has a clear +1/T dependence, the finely tuned 
singlets have a +1/T dependence at high temperatures 
that turns over to become a constant at zero temperature, 
resulting in a maximum at T ^ 2 A, where A is the sep- 
aration of the singlets. This temperature is significantly 
higher than the quenching of the linear susceptibility, Xi 
at T ^ .25A. By contrast, the Kramers doublet has a 
strong — dependence that should dominate at tem- 
peratures much smaller than the crystal field splitting, 

T<C Ayy. 

Of course, these are only single-ion calculations, and 
in URu2Si2 we only expect single-ion physics to hold 
for T > T* = 1{)K with (1/T ^ 1/(T - Ocw))^ 
making it difficult to distinguish the fine-tuned singlets 
from a true non-Kramers doublet. In the dilute limit, 
Ua;Thi_a;Ru2Si2 {x < .07), the single- ion physics ex- 
tends to much lower temperatures, down to lOK. Be- 
low this temperature, the Curie-like single- ion behavior 
is replaced by a critical logarithmic temperature depen- 
dence, - log T/Tk, where Tk - 10K[T5]. This physics 
has been attributed to two channel Kondo criticality, al- 
though magnetization studies show fine differences from 
the single impurity two-channel Kondo model [TO that 
may be due to the fact that two-channel impurities are 
never really in the idealized dilute limit. 

IV. DUALITY OF THE Fs DOUBLET 

As F5 is a magnetic non-Kramers doublet, it is pro- 
tected by a combination of both time-reversal and tetrag- 
onal crystal symmetry, (|±) = \x)^i\y))^ and it thus pos- 
sesses a dual magnetic and orbital character. This means 
that the the c-axis moment is magnetic (both dipolar and 
octupolar), while the basal plane moment is quadrupolar 
{Oxy,Ox'2_y2). As magnetic and quadrupolar moments 
result in linear and quadratic field splittings, respectively 

A(5„ B^) = Qf^pBc + ]^axvvBl, (19) 

they will appear in the linear and nonlinear susceptibil- 
ities respectively, as the nonlinear susceptibility is es- 
sentially the quadrupolar susceptibility. Thus the lin- 
ear c-axis susceptibility and the nonlinear basal plane 
susceptibility are predicted to have the same tempera- 
ture dependence - albeit with different coefficients. Tak- 
ing the above splitting, the free energy is given by F = 
— |A^, where r] = Xi/(^Mb)^ is found in c-axis field, 
F = —^Xi^c' In a basal plane field. 
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Figure 2: Linear susceptibilities for the three crystal field 
schemes with the c-axis and basal plane susceptibilities given 
by solid and dashed lines respectively. The basal-plane sus- 
ceptibility is multiplied by five, (a) Fs ground-state. Inset 
shows the anisotropic Curie behavior, (b) Finely-tuned sin- 
glets. Here the temperature has been rescaled by the gap, 
A = 36K between the two singlets. Inset shows the 1/T 
dependence of the magnetic susceptibility at temperatures 
larger than the gap A. (c) 5/^ Kramers doublet scenario. Inset 
shows Curie behavior of susceptibility in both the c-axis and 
basal plane. 



which implies that the basal plane nonlinear susceptibil- 
ity 

= 1.2 X 10"^ (T-f.u.)-\ (21) 
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Figure 3: The basal-plane nonlinear susceptibility, Xs" for 
the three scenarios, plotted versus the inverse temperature 
rescaled as T/A, where A = 36K is the gap for the singlet 
scenario. The Fs is given in solid red, the finely-tuned-singlet 
scenario in solid blue, and the 5/^ Kramers doublet scenario 
in dashed orange Note that only the singlet scenario has a 
strong non-monotonic temperature-dependence, with a max- 
imum around T = A. The F5 scenario has a strong, positive 
Curie-dependence, while the 5/^ is dominated by a negative 
-1/T^ term. 

will always share the temperature dependence of the lin- 
ear c-axis susceptibility, where we have used the crystal 
field parameters for the F5 doublet scenario (A). This 
shared temperature dependence is a result of the duality 
of the non-Kramers doublet and will hold not only in the 
single- ion limit, but in the critical dilute and dense lim- 
its. The Kramers doublet does not have this duality, as 
we can see from the different temperature dependences 
even in the single- ion case. The finely tuned singlets will 
only exhibit this duality down to the gap temperature 
(- 36K). 

V. DISCUSSION AND CONCLUSIONS 

Our treatment up to this point has assumed isolated 
single-ion behavior, determined by a nominal and inte- 
gral f- valence. In reality the 5f state is not a purely inte- 
gral valence state, but its magnetic properties do reflect 
a nominal integral valent configuration. The reasoning 
behind this is somewhat subtle and now deserves discus- 
sion. 

One of the important concepts that comes to our aid, 
is the distinction between the nominal integral valence 
that determines the magnetic properties and the micro- 
scopic valence measured by high energy spectroscopic 
probes pT]. To illustrate this point, consider the case of 
a nominal 5/^ configuration, undergoing valence fluctu- 
ations into a 5/^ singlet, 

5/' - 5/2 + e- 



In the isolated atomic limit, such a state is |5/^ : a), but 
once valence fluctuations are included, this state is then 
described by a Varma-Yafet wavefunctionp^ of the form 

|a*) = (l + ^afecL/a)|5f (22) 

k 

where the second term describes the virtual excitation of 
an f-electron into the conduction sea in a partial wave 
state with the same spin-orbit coupled angular momen- 
tum as the f-state. The important point here, is that 
while the valence is reduced below its nominal value, so 
that here nj- = 3 — (5, since valence fluctuations conserve 
total angular momentum, they do not renormalize the 
magnetic moment. For this reason, the magnetic moment 
of an isolated Anderson impurity is unaffected by signif- 
icant departures from integral valence [19 . The above 
arguments will of course also apply if the nominal con- 
figuration is a 5/^ non-Kramers doublet and Uf = 2 -\- S. 
The important point is that valence fluctuations do not 
renormalize the magnetic moment so that a single-ion 
treatment can be used to describe the susceptibility in 
the region that Curie behavior is observed. It also means 
that susceptibility measurements of any kind can only 
determine the nominal, not the microscopic valence of 
the ion. 

There is in fact a lot of circumstantial evidence that 
these kinds of arguments apply to URu2Si2. In particu- 
lar: 

• In the dilute limit, the magnetic susceptibility fol- 
lows a Curie Weiss form down to lOK^ giving us 
an estimate for the single-ion Kondo temperature 
Tk ~ lOK. In concentrated case, single-ion behav- 
ior continues down to about lOOK, which is still a 
small scale more characteristic of a system close to 
integral valence than a strongly mixed valent sys- 
tem. 

• STM measurements do indicate admixture between 
the f-electrons and the conduction electrons. In- 
deed, coherent tunneling into the f-states is ob- 
served on the Si layers of URu2Si2 , indicating 
that the renormalized f-states extend to the silicon 
orbitals[20il21j. 

Thus while there is mixed valence, single- ion magnetic 
physics is observed. 

One method to directly determine the microscopic va- 
lence of the U ions is from the branching ratios of two 
core- level d-states The quoted microscopic valence 
determined by these methods is n/ =2.6-2.8, which would 
suggest a nominal valence of 3 (5/^) with a Kramers 
doublet. Such measurements would at first sight seem 
to make a magnetic determination of the nominal va- 
lence unnecessary. However, recent theoretical work sug- 
gests that the microscopic valence determined by these 
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methods requires a more accurate many body treat- 
ment of the spin-orbit couphng in the f-states, which can 
give substantial corrections to the inferred microscopic 
valence [22]. 

In conclusion, we have proposed a bench-top exper- 
imental probe that can determine the nominal U va- 
lence configuration in URu2Si2, distinguishing between 
a 5f^ configuration that is naturally Ising-like and a 5f^ 
Kramers configuration that is fine-tuned close to the Ising 
limit. Using a single ion approximation valid at high 
temperatures, we have shown the basal plane non-linear 
susceptibility for a 5f^ non-Kramers or singlet-singlet 
ground-state is always positive, whereas the Kramers 
ground state turns negative at low temperatures. In the 
dilute limit, we have used the duality of the non-Kramers 
doublet to predict a temperature-independent ratio be- 
tween the basal plane non-linear susceptibility and the c- 
axis linear susceptibility; if found this will provide direct 
confirmation of the doublet. The nonlinear suscepti- 
bility of concentrated URu2Si2 was last measured twenty 
years ago [H [28]; measurements were only conducted for 
T < 25i^, well below the single-ion temperature regime. 
A set of new measusrements over an extended and higher 
temperature range in both the dilute and the concen- 
trated limits has the potential to add much insight to this 
problem. Resolution of the uranium single-ion ground 
state in URu2Si2 will provide crucial insight into the ori- 
gin of its hidden order, providing an important constraint 
on future microscopic theories. 
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(Flint, Coleman) and grant NSF 1066293 (Flint, Chan- 
dra, Coleman) while at the Aspen Center for Physics. 
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APPENDIX 

Here we explain the details of choosing a representative 
state for the 5f^ configuration. The strong spin-orbit cou- 
pling means that J = 9/2, and in a tetragonal field, the 
crystal field Hamiltonian contains five free parameters, 

H = BlO° + Biol + Biol + BtOl + B^Ot, (23) 

where are the Stevens operators [27]. There are three 
Tq doublets that mix = 9/2, 7/2 and 1/2 states and 
two Vj doublets that mix = 3/2 and 5/2, which be- 
tween them have eight free parameters, so the problem is 
grossly overconstrained. In order to find a set of param- 
eters, B that have the correct c-axis moment, high tem- 
perature anisotropy and van Vleck susceptibility, we per- 



formed a simulated annealing search of the space of 5's. 
We found that it is possible to find a state with nearly 
zero transverse moment, |(— |J+|+)| = .Ol/i^- This finely 
tuned state has a positive X3- (T) at low temperatures due 
to the lack of basal plane moment. However, this state 
also requires the existence of a low- lying doublet at llOK 
that is inconsistent with neutron scattering [1 and gives 
an entropy of greater than i?log4 at room temperature. 
And even in this case, X3"(^) is negative over some regime 
of temperature, so the sign of X3"(^) is still a good indi- 
cator of 5f^ behavior. 
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Figure 4: (Color online) Comparing possible 5f^ states: a 
finely tuned Ising state, .14| ± 9/2) + .39| ± 1/2) - .91| T 7/2) 
(purple, darker) and a more generic compromise state .05 1 =b 
9/2) + .42| ± 1/2) - .91| T 7/2) (orange, lighter). (A) The 
linear susceptibilities both look similar at high temperatures, 
though the absence of a basal plane moment can be seen at 
low temperatures. (Inset) Crystal field scheme for the finely- 
tuned state. (B) The finely-tuned state has a low-lying dou- 
blet at llOK that makes the low temperature entropy much 
larger. (C) The nonlinear susceptibility looks the same above 
lOOK, but the lack of basal plane moment gives a positive 
contribution to the finely-tuned state that resembles the Fs 
upturn. However, the negative X3"(^) above lOOK clearly in- 
dicates the 5f^ nature. 



